Dynamical spin-flip susceptibility for a strongly interacting ultracold Fermi gas 
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The Stoner model predicts that a two-component Fermi gas at increasing repulsive interactions 
undergoes a ferromagnetic transition. Using the random-phase approximation we study the dynami- 
cal properties of the interacting Fermi gas. For an atomic Fermi gas under harmonic confinement we 
show that the transverse (spin-flip) dynamical susceptibility displays a clear signature of the ferro- 
magnetic phase in a magnon peak emerging from the Stoner particle-hole continuum. The dynamical 
spin susceptibilities could be experimentally explored via spin-dependent Bragg spectroscopy. 
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Introduction Experimental advances of trapping and 
cooling ultracold Fermi gases allow the exploration of 
strongly interacting quantum degenerate two-component 
gases, realized by populating two hyperfine atomic lev- 
els. The strength and the sign of the intercomponent in- 
teractions can be tuned by the mechanism of Feschbach 
resonances. The case of attractive Fermi-Fermi interac- 
tions has allowed the exploration of the crossover from 
BCS-like pairing to Bose-Einstein condensation (BEC) 
of tightly bound pairs On the repulsive side of the 
Feschbach resonance, the metastable branch correspond- 
ing to a two-component Fermi gas with repulsive inter- 
actions has received more recent experimental attention 
0, [3 as a possible realization of the Stoner model of itin- 
erant ferromagnetism. The mean-field Stoner model [4[ 
predicts that a repulsive Fermi gas undergoes a transition 
to a ferromagnetic state for sufficiently large interaction 
strength kpa > tt/2, where a is the s-wave scattering 
length and kp — (S^n) 1 ^ 3 is the Fermi wavevector for 
a uniform gas with total density n. Beyond mean field 
corrections, the inclusion of quadratic fluctuations de- 
creases this value to kpa » 1.05 0], and Quantum 
Monte Carlo simulations predict the transition to occur 
at kp a » 0.82 0, Hj. While these calculations assume a 
uniform ferromagnetic phase, the possibility of inhomo- 
geneous intermediate phases has been considered Q. 

In this work we analyze the dynamical properties of 
a two-component Fermi gas with repulsive interactions. 
The approach to the ferromagnetic transition strongly 
affects the dynamical properties of the system, display- 
ing a softening of the out-of phase spin-dipole modes [10( 
and enhancement of the spin drag coefficient [11]. A dy- 
namical instability of a spin spiral excitation has also 
been predicted [12[. Above the ferromagnetic transition, 
magnon collective modes are predicted to propagate un- 
damped at zero temperature |13| . 

Spectroscopy is a powerful tool to address many-body 
systems. The study of the dynamic structure factor for 



density and spin has been suggested as a probe of pair- 
ing and superfluidity of a Fermi gas [l4, 15|. and Bragg 
scattering experiments for a strongly interacting Fermi 
gas in the BCS-BEC crossover have been performed [lij ]. 

We focus here on the spin-susceptibility spectra as a 
tool to characterize the ferromagnetic phase of a two com- 
ponent Fermi gas. Such spectra could be accessed exper- 
imentally via spin-dependent Bragg spectroscopy. 

Dynamical susceptibilities in the random phase approx- 
imation We consider a two-component Fermi gas | f), 
\ \), at zero temperature, with repulsive intercomponent 
interactions described by a contact potential of strength 
g = Anh 2 a/m in terms of the scattering length a > 0, and 
the atomic mass m. The gas is subjected to an external 
confining potential V ex t(r). The model Hamiltonian of 
the system reads 



H 



2m 



V 2 + V ext {r) 



Vv(r) 



; a a are the Pauli matrices with a = 0,x,y,z. 



+g J drV4(r)^(r)^(r)V> t (r), (1) 

where the fermionic field operators V'j.(r), - 0t( r ) satisfy 
the usual anticommutation relations {Vv(r), W a i ( r )} = 
S a . a 'S(r — r'). The corresponding spin densities are 
CT Q (r) = v P t (r)cr Q ^'(r) in terms of the spinor \&(r) = 
( V>t(r) 

For the transverse (i.e. x, y) components it is useful to 
introduce the combinations a± = {a x ± io~ y )/2 because 
they yield the spin-flip amplitudes 

a + (r) = V4(r)^(r); a.(r) = $(r)^(r). (2) 

In linear response theory, corresponding to the regime 
of weak external drive, the dynamical spin susceptibilities 
are given by the retarded Green's functions 

Xa p(r,r',t,t') = -i6(t-t')([o- a (r,t),*p(r',t')}), (3) 
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where the quantum averages (...) are performed on the 
unperturbed state of the system. Such dynamical cor- 
relations functions embed the effect of interactions, and 
require a many-body calculation. We consider here their 
expression as obtained by the Random-Phase approxi- 
mation (RPA) [I?], EH , corresponding to time-dependent 
Hartree-Fock approximation. This analytical approxima- 
tion is known to satisfy sum rules and well describe the 
weak-coupling regime. It also accounts for the ferromag- 
netic transition at the mean-field (Stoner) level, hence, 
although not quantitatively correct it is expected to qual- 
itatively describe the dynamical properties of the inter- 
acting Fermi gas. 

According to the RPA, for a uniform Fermi gas the 
Fourier transform of spin-flip susceptibility with respect 
to both spatial and time relative variables is given by 



x+-( w >o) = 



IF 



1 + 5 IL, 



(4) 



A similar result holds for the total density xoo an d 
longitudinal-spin \zz responses, 



Xoo(w,q) 



n ++ + n__ + 2.gn ++ iF 



Xzz{u-,q) = 



1 - ,g 2 n ++ n_ 



n ++ + n__ -2 5 n ++ n_ 



(5) 



i-. g 2 n ++ n__ 

where the Lindhard functions II CTCT '(u;,q) are given by 
IW (w, q) - - ^ Z—JP -7f—„ 



(6) 



(7) 



with er, a' = ± indicating the spin polarization along the 
axis of spontaneous magnetization, f a (Q = (e"" + 
(3 = 1/ksT and = — \i + gn- a . We have used 
the symbol ± to indicate the two magnetization compo- 
nents since for long transverse decoherence times and in 
absence of population imbalance the direction of magne- 
tization of the ferromagnet is not necessarily the quanti- 
zation direction of the total angular momentum F of the 
atoms used to define the two-component Fermi gas. 

Below the ferromagnetic transition, ie in absence of 

population imbalance we have 11^ = II ++ = II 

1 1 . the well-known Lindhard function for noninteracting 
fermions [ljj], and Xoo(zz) = 211/(1 Tffll). For repulsive 
interactions, characterized by g > 0, these expressions 
predict a well-defined propagating collective mode for the 
density channel (the zero sound) and a damped mode in 
the spin channel. 

The validity of the RPA can be inferred by compari- 
son with the Landau theory of Fermi liquids. The latter 
allows to calculate the density (and spin) susceptibility 
Xoo(zz) a t small energy and long-wavelength below the 
ferromagnetic transition. The obtained susceptibilities 
[20l | have the same structure as Eqs.© and ©, with 




FIG. 1. (Color online) Imaginary part of the spin-flip sus- 
ceptibility Inixn (q, w) for a uniform Fermi gas in the (w, q) 

plane (q in units of %f, w in units of hkp/m) for the fully 
polarized case n+ = n. The magnon dispersion is indicated 
by the solid white line. 



the interaction strength g replaced by suitable Landau 
parameters Fq^. 

Spin-flip susceptibilities for a uniform Fermi gas above 
the ferromagnetic transition 

Above the ferromagnetic transition a window opens in 
the (w, q) plane at small frequencies and wavevectors, 
below the Stoner particle-hole continuum, allowing for 
the propagation of spin- flip collective modes (magnons), 
corresponding to the Goldstone mode of the ferromag- 
netic phase. The magnons have a quadratic disper- 
sion n q which is obtained from the poles of Eq.@, ie 
1 + gll_| (rt q , q) = 0, yielding at second order in q 



5A 2m 



(n+fc 



-kr 



(8) 

with A = g(n+ — ri—), recovering [13j in the fully polar- 
ized case, and kp = (6ir 2 n±) 1 / 3 . 

The contribution of the particle-hole continuum is ob- 
tained from the study of the Lindhard function II_| , for 

which we obtain an analytic expression. For simplicity 
of notations we rescale |q| in units of kp = (k + 



and set v 



muj/hk F , 5k = (kp 



k F )/2k F 



mg(n + — n )/hkp. The Lindhard function at zero tem- 
perature is nonzero only in a restricted part of the (y, q) 
plane, delimited by the curves u^' 2 ' = q 2 /2+d±q(l+Sk), 

= —q 2 /2 + d =F ?(1 — Sk). In detail, by setting 
A± = (1 ± 6k), B± = ((v - d)/q =p <?/2), in the rescaled 
units rakp/2h 2 for II_| the imaginary part reads: 



QTI , _ = 



-1 



A 2 



B\ - A 2 _ 



B 2 



in (a) 
in (b) 
in (c) 



(9) 
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where the regions in the (v, q) plane are defined by the 
following inequalities: (a) max{z^_ 2 \ v^} < v < or 
(b) max.{i/^\ v+^} < v < i/^; and (c) 

(2) ,,(2)^ 



(2) , , (1) 



part we have 



v / <v< min{i/\ , v_ }. Correspondingly, for the real 



2ir 2 q 



1 



[A+B+ - AS. 
A + + B + 



Si) In 



A+ - B. 
A- + B_ 



(10) 



A_ -B- 



The above results for the Lindhard function allow to 
obtain the final result for the imaginary part of the spin- 
flip susceptibility 



SIl + -(u,q) 



(i + 5 5Rn+_( Wj(Z ))2 + ff 2 3 2n + _ 

+A(q)S(u - fi g ), (11) 



this includes both the Stoner particle-hole continuum 
Sx+™(w, q) and the magnon contribution with weight 

A{q) = (g 2 \dll + -/d U \ at )- 1 - 

The RPA expression for the spin- flip susceptibility sat- 
isfies the sum rule [2l| 

- - ^ duZ X+ -(uj,q) = ( n +~ n ~) . (12) 

This provides an alternative way to fix the weight A{q) 
of the collective magnon mode with respect to the con- 
tribution of particle-hole excitations, according to the re- 
lation A(q) = (n+ - n J)/h + (1/tt) C x duSSx^i", <?)• 
At increasing q values the magnon strength decreases till 
it vanishes when the magnon mode enters the particle- 
hole continuum. Nevertheless, as it is shown in FigfTJ 
quite remarkably some trace of the magnon mode is left 
as a pronounced maximum in the particle-hole contin- 
uum. This feature could be used to infer the presence of 
the magnon mode even for wavevectors larger than those 
associated to the Stoner gap. 

Spin-flip response for a trapped gas 

We consider now the effect of inhomogeneity in the 
Fermi gas induced by the presence of the external har- 



monic confinement V ex t(r) 



taken for sim- 



plicity isotropic. The equilibrium density profiles n+(r), 
n_ (r) of the two magnetization components are found by 
minimizing the Thomas- Fermi energy functional [22I [23j | 



E[n a (r)] 




1 



cm CT (r) 5/3 + -mttt r r 2 n a (r) 



5 v ' 2 
ff n+(r)n_(r)} 



(13) 



where a = (67r 2 ) 2 / 3 fi, 2 /2m. The chemical potentials ii a 
are related to particle numbers through the normaliza- 
tion conditions J drn CT (r) = N a . Note that for simplicity 
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FIG. 2. (Color online) Main panel: Spatial density pro- 
files n+(r) (top) and n_(r) (bottom), in units of the cor- 
responding noninteracting densities at the trap center n^ ' (0) 
as a function of the radial distance r in units of auo for 
dimensionless interaction strength A = k%(0)a = 2.5, with 
k%(0) = (67r 2 nS 0) (0)) 1/3 , and for a particle number N = 10 7 . 
Insets: the gray area indicates the Stoner particle-hole contin- 
uum in the (to, q) plane as seen at various regions in the trap 
(center, middle, and sides from left to right) in a local-density 
picture, and the red circle indicates the region where magnon 
excitations are undamped at zero temperature. 



the calculations have been performed at = fi-, cor- 
responding to fix the total number of particles N. An 
analysis beyond Thomas- Fermi approximation [24j would 
provide corrections at the tails of the density profiles but 
not alter significantly the resulting profiles. 

The dynamical spin-flip response of the trapped sys- 
tem is obtained from the one of the homogeneous system 
through the local-density approximation (LDA) as an av- 
erage over the various trap regions, 



Iv tr rfr g s x+-(^,q,r) 

fv tr dr 



(14) 



where Sxh (w, q, r) is the homogeneous result (fTTj) 

which depends on r through n+(r),n_(r), and Vt r in- 
dicates the spatial volume where the particle densities 
are non vanishing. The local-density approximation is 
valid if the typical excitation frequency and wavevectors 
are larger than harmonic-oscillator frequency f2t r and the 
inverse size of the cloud. Figure [5] illustrates the particle 
density profiles and the Stoner particle-hole continuum 
corresponding to different trap regions: while in the trap 
center the Stoner gap is maximal and magnon propaga- 
tion is allowed, at the trap sides the gap closes and only 
the particle-hole continuum contributes to the dynamical 
response. 

We illustrate in Figj3]the total spin-flip response as a 
function of the frequency, for values of system parameters 
accessible in current experiments. Even in the presence 
of the external confinement we find that at small trans- 
ferred wavevectors the dynamical response is dominated 
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FIG. 3. (Color online) Imaginary part of the spin-flip sus- 
ceptibility — glmx+-(q, w) for a harmonically trapped Fermi 
ff£LS £LS El function of the rescaled frequency Q = muj/hk ( p(0) 2 
at various values of the transferred momentum Q = q/k%(0) 
for dimensionless interaction strength A = 2.5. The hetched 
areas indicate the magnon collective excitation contribution 
to the spectrum, the solid lines indicate the contribution of 
the Stoner continuum which displays a maximum when the 
magnon merges onto it. 



by the magnon contribution, emerging from the particle- 
hole continuum. The latter acquires more importance 
at larger values of wavevectors. We expect this picture 
to hold at finite temperature, provided that it is smaller 
than the Stoner excitation gap A. 

Conclusions In this work we have considered the dy- 
namical spin-flip response of a repulsive Fermi gas above 
the ferromagnetic transition, which is characterized by a 
magnon collective excitation mode. Within the random- 
phase approximation we have provided an analytic ex- 
pression for the the spin-flip response of the homogeneous 
two-component Fermi gas. For the experimentally rele- 
vant situation of an inhomogeneous gas, we have shown 
that even in the presence of an external confinement 
the dynamical response displays features of the magnon 
mode in a pronounced peak, thus displaying a clear in- 
dication of the ferromagnetic phase. These predictions 
are experimentally accessible via spin-dependent Bragg 
spectroscopy, which could be realized, similarly to the 
usual Bragg spectroscopy, by a two-photon process in- 
ducing transitions from | fermions to | 1) fermions, 
transferring at the same time momentum hq and energy 
hui to the fluid. 

Refinements of our model include the development 
of a fully quantum description for the spectrum of the 
trapped interacting gas both at RPA level as done for 
the paramagnetic phase in 



tion would require to include in the dynamical descrip- 
tion the effects of atom losses, and the presence of bound 
states as in [27 1. 
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